Preprint typeset in JHEP style - HYPER VERSION 



Bispectrum signatures of a modified vacuum in single 
field inflation with a small speed of sound 



P. Daniel Meerburg 1 , Jan Pieter van der Schaar 2 and Mark G. Jackson 3 

1 Astronomical Institute "Anton Pannekoek", University of Amsterdam, 
Science Park 904, 1098XH Amsterdam, The Netherlands 

2 Korteweg-de Vries Institute for Mathematics, University of Amsterdam, 
Plantage Muidergracht 24, 1018 TV Amsterdam, The Netherlands 

3 Lorentz Institute for Theoretical Physics, 2333 CA Leiden, The Netherlands 



Abstract: Deviations from the Bunch-Davies vacuum during an inflationary period can 
leave a testable imprint on the higher-order correlations of the CMB and large scale 
structures in the Universe. The effect is particularly pronounced if the statistical non- 
Gaussianity is inherently large, such as in models of inflation with a small speed of sound, 
e.g. DBI. First reviewing the motivations for a modified vacuum, we calculate the non- 
Gaussianity for a general action with a small speed of sound. The shape of its bispectrum 
is found to most resemble the 'orthogonal' or 'local' templates depending on the phase of 
the Bogolyubov parameter. In particular, for DBI models of inflation the bispectrum can 
have a profound 'local' template feature, in contrast to previous results. Determining the 
projection into the observational templates allows us to derive constraints on the absolute 
value of the Bogolyubov parameter. In the small sound speed limit, the derived constraints 
are generally stronger than the constraint obtainable from the power spectrum. The bound 
on the absolute value of the Bogolyubov parameter ranges from the 10 -6 to the 10 -3 level 
for H/A c = 10 -3 , depending on the specific details of the model, the sound speed and the 
phase of the Bogolyubov parameter. 
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1. Introduction 

Recent years have witnessed a steady increase in the amount of precision data from the 
cosmic microwave background and large scale structure [1]. This opportunity has moti- 
vated the study of both the power spectrum (or two-point correlation) of the primordial 
density perturbations, essentially probing the free theory of the inflationary Lagrangian, 
as well as higher-order statistics, probing the interacting theory. The leading statistical 
probe of interactions is the bispectrum (or three-point correlation, indicating non-Gaussian 
statistics) and there is currently a large collection of work on the non-Gaussian predictions 
for specific models of inflation [2-6]. Interestingly, any detection of non-Gaussianity will 
immediately rule out all models of single-field slow-roll inflation. While this would be a 
monumental discovery in itself, it would also be the starting point for further constraining 
the inflationary Lagrangian. Indeed, if the non-Gaussian amplitude is high enough to be 
observable it opens up a new, and very distinctive, window on the details of inflation [7] . 

In a previous article by two of us (PDM, JPvdS) [8] we computed the primordial bis- 
pectrum for a single-field slow-roll inflationary model, with the assumption that the initial 
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state is not the Bunch-Davies (BD) vacuum [9]. This assumption is well-motivated because 
inflation is by construction an effective theory and therefore there is no 'knowledge' of any- 
thing beyond a specified cutoff scale in energy In expanding backgrounds this energy-cutoff 
is tantamount to a (possibly momentum-dependent) time-cutoff, allowing one to instead 
parameterize ignorance as an initial state from which to evolve the dynamics. As with 
all effective theories, the details underlying the construction of such a state are irrelevant; 
the only concern is to use this state to compute the consequences on the bispectrum, and 
if possible compare theoretical predictions with the data to place constraints on such a 
deviation. 

As first reported in [10], our results in [8] confirmed that small initial state modifi- 
cations can induce surprisingly large corrections to the primordial bispectrum, especially 
if one includes higher-derivative operators (see also [11]). The intuitive explanation is 
essentially that any deviation from BD introduces particles on sub-horizon scales at the 
initial stages of inflation, and these particles then have ample time to produce a large non- 
Gaussian signal in the presence of (irrelevant) interactions. Unfortunately, this yielded 
only a rather weak bound of order |/3| < 1CT 2 on the deviations from the BD state in 
terms of the Bogolyubov /3-parameter. This is because although the magnitude of such 
induced non-Gaussianity was large, its shape in co-moving momentum space was found 
to substantially depart from standardized observational templates, providing poor overlap 
with existing constraints. This motivates the need for better templates and possible sys- 
tematic analysis of oscillatory features in the bispectrum, for which we hope to report some 
progress in the future. 

In this article we generalize the class of models for which initial state modifications are 
considered. We opine that initial state modifications are generally natural for these models. 
In particular, we will argue that for both branches in DBI inflation, UV [12] and IR [13,14], 
small deviations can be expected or at least are worth considering when computing the full 
bispectrum. Our starting point will be the most general single-field action, for which the 
bispectrum has been previously computed in [15, 16]. Specific single-field models (slow-roll, 
DBI, kinetic) are all limiting cases of this action, and we will be interested in the limit 
where the speed of sound c s (that is, the speed at which inflaton perturbations evolve) is 
significantly reduced relative to c, the speed at which metric perturbations evolve. Whereas 
bispectrum components in standard slow-roll models are typically proportional to e and 
r\ [17], in the limit of small c s one finds new contributions to the bispectrum that are 
proportional to l/c 2 s . Assuming the BD- vacuum these potentially large corrections are 
typically of the equilateral type and have been previously constrained [18, 19]. When 
one allows for (initially Gaussian) vacuum state modifications, there are yet-additional 
components to the bispectrum which are also expected to be enhanced in the small speed 
of sound limit. As a consequence, we will be able to put relatively strong constraints on 
\(3\. Importantly, we will also study the effect of an arbitrary phase in the Bogolyubov 
parameter, which has so far been neglected. Somewhat surprisingly, we will find that 
specific choices for the phase can have dramatic effects on the shape of the leading non- 
Gaussian contribution. As a consequence the phase of the complex Bogolyubov parameter 
(5 strongly affects the observational constraints on the absolute value \(5\ that can be derived. 
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An important restriction in our analysis is that we will have to assume approximate scale- 
invariance throughout and will not discuss (big or small) departures from scale invariance 
and how these can be observationally constrained [20,21]. 

This article is organized as follows. In §2 we will briefly summarize the small-c s models 
of inflation and explain how initial state modifications are quite natural in this context. 
In §3 we begin from the most general single-field inflationary action and then compute the 
corrections to the bispectrum. In §4 we identify the leading-order and sub-leading com- 
ponents and then constrain modifications to the initial state. Constraints can only be set 
using already measured bispectra by comparing the theoretical components to observed 
templates. In order to achieve this, we compute the leakage via projection of these com- 
ponents into measured components, for which we make use of a formalism first proposed 
in [22] and optimized in [23]. In §5 we use the two derived constraints on different /nl 
to derive constraints on DBI models with three free parameters. In §6 we summarize and 
conclude. 

2. A case for vacuum state modifications in small sound speed models 

It has now become clear that vacuum state modifications cannot be excluded on purely 
theoretical grounds, even though they are constrained by bounds on backreaction [24, 25] 
and observations of the primordial power spectrum [26-30]. The theoretical challenge is to 
relate potential vacuum state modifications to ultraviolet physics at the string- or Planck- 
scale. In terms of an effective field theory description of inflation such a relation would 
appear natural, since the details of the vacuum (or equivalently initial) state of inflation 
usually depends on physics far beyond the cutoff in the effective field theory description. 
This raises the hope of being able to probe string- or Planck-scale physics through the 
signatures of possible modifications to the vacuum state. 

Unfortunately, so far the available proposals have been phenomenological in nature and 
can be separated into two classes. In the original work of Danielsson [31] scale-invariance 
is preserved by proposing that the initial time rjo (beyond which an effective field theory 
theory breaks down) is identified as the time where the physical momentum equals some 
high energy cutoff, defining a New Physics Hypersurface (NPH). This immediately implies 
that the initial "cutoff" time will be different for different comoving momenta, i.e. rjo will 
be a function of k allowing for the preservation of scale-invariance. A specific, but rather 
ad-hoc, prediction for the Bogolyubov parameter is then obtained by assuming a local 
empty state at the initial time rjo(k) [31,32]. 

In a different approach known as Boundary Effective Field Theory (BEFT) [33], scale- 
invariance is explicitly broken by identifying an initial time surface rjo where the usual 
BD boundary conditions are corrected by including leading corrections. As one would 
expect from effective field theory, the corrections to the BD state grow as a function of 
comoving momentum k, with the largest comoving momenta corresponding to physical 
momenta closest to the cutoff scale on the initial time surface. The BEFT approach has 
the advantage of allowing one to calculate corrections to the BD state systematically from 
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first principles. However, the breaking of scale-invariance and the dependence upon the 
initial time r/o are serious drawbacks. 

In both scenarios the magnitude of the (leading) corrections, parameterized by the 
absolute value of the Bogolyubov parameter, is proportional to the inflationary Hubble 
parameter H divided by a physical cutoff scale A c . The appearance of the Hubble scale 
is natural for various reasons and can in particular be related to the fact that an initial 
vacuum state can only be defined in the adiabatic regime, i.e. for physical momenta much 
bigger than the Hubble scale. This immediately suggests that any modification should be 
proportional to H/A c , ensuring that the breakdown of a small (perturbative) initial state 
correction coincides with the breakdown of adiabatic behavior. 

This brings us to our first general observation regarding inflationary models with a 
small speed of sound. As a direct consequence of the small speed of sound the adiabatic 
regime breaks down long before the Hubble scale is reached. Instead, the relevant scale 
in this case is H* = H/c s 3> H, whose inverse corresponds to the radius of the sound 
horizon. Consequently, a vacuum state modification is expected to be proportional to H* 
in any small speed of sound scenario. Therefore irrespective of a particular motivation or 
method to describe vacuum state modifications in models with a small speed of sound, 
its magnitude is anticipated to be enhanced by a factor of l/c s as compared to standard 
slow-roll inflationary models (for the same Hubble parameter). Taking into account the 
constraint on c s from bispectrum data, which roughly corresponds to c s > 10~ 2 [19], 
this could add two orders of magnitude to the observable range of H/A c as compared to 
standard slow-roll inflation. 

The identification of the sound horizon as the relevant scale, signaling the breakdown 
of the adiabatic approximation, also has a more fundamental consequence for the initial 
vacuum state. Generically, the speed of sound slowly evolves during inflation and the 
models can be separated in two classes: ones where the speed of sound decreases (as in UV 
DBI models) or models where the speed of sound increases (as in IR DBI models) towards 
the end of inflation. In the second case, as was first pointed out in [34] , when e > 1 — s 
(where e is the usual slow-roll parameter and s = c s /H c s measures the rate of change in 
the speed of sound) the fluctuation modes evolve from superhorizon to subhorizon, which 
excludes the possibility of an adiabatic initial state such as the Bunch-Davies state. Of 
course, this also rules out any phenomenological approaches to vacuum state modifications 
that require an (approximately) adiabatic regime, such as the NPH and BEFT proposals. 
Instead the (different) initial vacuum state in this case should depend crucially on the 
physics before the inflationary phase. In IR models of DBI inflation this would correspond 
to the region near the tip of the throat where (strongly coupled) stringy physics is expected 
to become important. 

In the opposite scenario, where the sound speed is decreasing towards the end of 
inflation, there is also an additional reason why the vacuum could be different from the 
standard BD vacuum [34]. In these models the speed of sound at the time of horizon 
crossing of the scales observable today should be large enough (roughly c s > 10 -2 ) to avoid 
violating current bounds on non-Gaussianity. Depending on the values of the "slow-roll" 
parameter e and the sound speed rate of change s, this requirement can imply that there 
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Figure 1: DBI inflation; IR and UV branches resulting in different evolution for c s , leading to 
different motivations for changes in the Initial Conditions (I.C.) 

are not enough e-folds left for the largest observable scales to allow for an approximately 
adiabatic vacuum choice. In other words, some modes in this case are not in the adiabatic 
regime when inflation began and are probing, just as in the previous case, the phase before 
inflation started. In the specific example of UV models of DBI inflation (whose simplest 
versions are in fact already ruled out), these modes would be probing the bulk region 
outside of the throat, which is governed by stringy physics (see figure 1). 

Finally, considering in particular DBI models of inflation, it is important to realize 
that the DBI action is obtained as the effective theory of open strings after heavy closed 
strings have been integrated out. As was pointed out in [35] integrating out heavy (closed 
string) states during an inflationary stage could give rise to modified initial conditions in 
the effective theory of the light (open string) degrees of freedom. 

The arguments above provide a relatively strong case for considering departures from 
the standard BD vacuum in models with a small speed of sound, in addition to the generic 
effective field theory motivation. It should therefore be worthwhile to look for bispec- 
trum signatures of modified initial states in the limit of a small sound speed, especially as 
compared to canonical slow-roll inflation. In our analysis we will consider initial vacuum 
state modifications parameterized by an arbitrary complex Bogolyubov parameter, so no 
particular model for vacuum state modifications is assumed. As first reported in [10], and 
corroborated in [8, 11], the bispectrum is significantly enhanced due to the presence of par- 
ticles in the modified initial vacuum state. The important requirement of scale-invariance 
in order to do our analysis, in addition to the need for a (comoving dependent) cutoff 
scale to regulate the integrals after a change in the initial vacuum state, implies that the 
enhancement factor is expected to be proportional to a power of the ratio This is 
explains why the bispectrum can give rise to relatively strong constraints as compared to 
the power spectrum. Even though the appearance of H* = H/ c s implies that the enhance- 
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ment is reduced (for fixed H) in the small sound speed limit, the specific powers of the 
enhancement factor that will appear, together with the expected overall enhancement of 
the non-Gaussian signal in the small sound speed limit, will result in a bispectrum whose 
magnitude is significantly larger as compared to the c s = 1 case. 

3. Corrections to the bispectrum 

3.1 General single-field models of inflation 

We will be working in the context of a general (minimally coupled) 4d action 

S = I J d 4 x^[M 2 R - 2P(X, eft)], (3.1) 

where (ft is the inflaton field and X = g^ u d^(ftd u (ft. From this action we can obtain expres- 
sions for the usual slow- roll parameters, as well as a few parameters relevant for the three- 
point function first introduced in [15, 16] which are proportional to higher-order derivatives 
of P(X,<f>): 

XP, X 
M 2 H 2 ' 

P X = MpPe 

P x + 2XP XX £ ' 

XP X +2X 2 P XX , 

X 2 P XX + 2 -X^P^ xxx = l - (x^ - . (3.2) 

In a recent paper by Senatore et al. [19] the description of single-field inflationary fluc- 
tuations was generalized further by focusing on the (broken) time-translation symmetry 
and writing down the most general effective Lagrangian of the corresponding Goldstone 
boson degree of freedom. In their description contributions to the three-point function are 
organized differently and known results derived in the standard formalism are only repro- 
duced after appropriately translating the relevant variables. In particular they introduce 
a parameter C3 relevant for the (shape of the) three-point function, which is related to the 
above set of parameters as 



As pointed out in [19], whereas the generic expectation for the shape of the three-point 
function in single-field models with a small speed of sound is equilateral, i.e. a maximal 
signal in equilateral triangles, at special values of C3 (around —5) the non-Gaussian signal 
peaks in enfolded (or collinear) triangles instead. In fact, the orthogonal template put 
forward in [19] to analyze the CMB data is closely related to an earlier proposal in [8] 
designed with initial state modifications in mind, which also tends to peak in enfolded 
triangles. It would be interesting to find a specific Lagrangian function P(X, (ft) that 
can lead to these special values of C3, which might be impossible. Subsequently we will 



£ = 
A = 



2A 
m ~ c 2 s ) 



(3.3) 



- 6 - 



base our three-point function calculations on the more "traditional" general single-field 
action described by a function P(X, </>), which can always be mapped to the corresponding 
Goldstone boson action using (3.3). 

A specific example of much interest is the string-inspired model of DBI inflation [12, 13] , 



P(x,</>) = _/(0)-yi- 2A7(0) + f{4>r 1 - v{4>). 

In this case, 

H 2 M% e 

Substituting this into (3.3) one finds the corresponding C3 = |(1 — c? s ), which equals 3/2 
in the small sound speed limit. This will cause the leading contribution in the three- 
point function to vanish for DBI models, which will be of importance when comparing the 
theoretical non-Gaussian amplitudes and shapes to the observational constraints. 

In order to use the action given by (3.1) to calculate our quantity of interest, the 
three-point function, we need to determine the corresponding interaction Hamiltonian for 
the gauge-invariant perturbation variable. After first perturbing in the inflaton field <j> = 
4>q + 5(j), we then convert this quantity to the gauge-invariant fluctuation £ = — (H/<po)6<j). 
Expanding the action to third order in these perturbations, the effective interaction Hamil- 
tonian associated with the cubic action is [14] 

Hj(f,) = -j d 3 x -a (E(l - 1) + 2A) ^ + ^(e - 3 + ?,c 2 s )Ul 
2 

+ °Ll{e-2s + l-c 2 s )UdQ x ) 2 -2^aUd( x )(d X ) , (3.4) 

where over dots represent derivatives with respect to conformal time. This interaction 
Hamiltonian will serve as our starting point for the perturbative calculation of the three- 
point function in a modified initial state. 

3.2 The three-point function 

We are interested in calculating the three-point correlation function {CfciCfc 2 Cfc 3 ): where Cfc 
is the gauge-invariant primordial density perturbation of comoving momentum k. To first 
order in the interaction Hamiltonian Hj this correlation can be written as an integral over 
the free-field correlator (Ck 1 Ck 2 (k 3 H[(r])) , 

(miCkAvKkMCkMmv)) = (Hm)\(kMCkM(kM\Hm)) (3.5) 



where the time-evolved initial state is given by 

mrj))=Te- lI "o HlW)d71 'm m )). 
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The first term on the right in eq. (3.5) will be zero if the initial state \ip(rjo)) is Gaussian, 
as we will assume. The interaction Hamiltonian will encode all higher-order contributions 
in the effective action, so it will be sufficient to consider interactions to lowest order in Hj. 
Recall that if the initial state is assumed to be BD at the onset of inflation, the 3-point 
correlator at tree level is usually computed from rjQ — > — oo. In our case we assume there 
is some part of the physics of which we are ignorant, ending up in the initial state at the 
start of inflation. Consequently we integrate from some initial time r/o representing the 
initial time corresponding to the cutoff scale in physical momenta p = k/a. To stress this 
again, 770 will therefore be fc-dependent and lead to a scale-invariant spectrum, which we 
require for our analysis. The presence of such a cutoff must be explicit because initial state 
modifications would otherwise have infinite time to interact, and the 3-point correlator 
would then be infinitely large - and hence excluded by observation! In our case it is the 
time at which we chose to set the initial state to be Gaussian-distributed parameterized 
by an a priori arbitrary complex Bogolyubov parameter. Relating the initial time to the 
moment the physical momentum of a mode reaches the cutoff scale implies that all modes 
have equal time to evolve until they exit the horizon, which allows for the preservation of 
scale-invariance in the bispectrum. 

One can expand the second line of the three-point correlator (3.5) into a product of 
Wightman functions G^. These are defined as 



The Wightman functions can be found by solving the classical equations of motion of the 
inflaton field action eq. (3.1) minimally coupled to gravity, 



The quantized form of the dimensionless curvature perturbation is given by 

Ck(v) = u k (r])a k + u*_ k (vW- k , 

where a k and at are the annihilation and creation operators respectively and the mode- 
functions Uk{rj) are dimensionless by introducing the appropriate powers of the reduced 
Planck mass M p . Using this form of Ck(v) one obtains the Wightman function 




(3.6) 




iH 



ike- 




H 2 



We compute the late-time cases of interest to us, 



G£(0,7,) 



H 2 



(1 



ikc s r])e 



4M 2 ec s k 3 



and 




(0,r?) 



H c s 



s ikCeV 



(3.7) 



4M 2 ek a{r{) 
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with a (77) = —1/rjH during inflation in the usual assumption H ~ 0. To explicitly compute 
the correlation (3.5), a useful identity is 



{CkMCkMCkM) 



-21Ze 



idrj'i^lCkMCkMCkA^HimiP^o)) 



, (3i 



which removes the commutator and assures that we are dealing with a true observable. 
Next we can express {Ck 1 Ck 2 Ck 3 ) in terms of these Wightman functions. For illustrational 
purposes let us just consider the first term on the right in eq. (3.4) in the absence of initial 
state modifications. The three-point function at late times 77 — > is then 



(Cfci(0)Cfe,(0)Cfc3(0))BD = -12(27r) 3 «5( 3 )(^^ 



ii-4i 



Tie 



1 



i / dr,' adjGlMdtfGlMdrfGl&r,') 



The factor 12 comes from 2 times the real part and the 6 possible permutations, and we 
have added the subscript to indicate that this is the result for the Bunch-Davies vacuum. 
From here on we will omit the 77 = dependence of the correlation function, since we are 
only computing it at late times. Using the Wightman functions (3.7) we obtain 

,3 r r0 dr] > 



S 1 



1 



+ 2A 



8M p 6 e 3 



Tie 



1 



ik t c a r)' 



(3.9) 



where kt = k\ + + £3. In the limit 770 
becomes 







dr]' ri' 2 e iktCsV ' 



m « 2 hk 2 k 3 
00, the (analytically continued) integral 

H 2 2 



hk 2 k 3 

Substituting this into eq. (3.9) we find 

(CkMk^BD = -3(2vr) 3 <5( 3 )(^ 



k}kik 2 k 3 zc 3 



S I 



+ 2A 



H 2 



1 



8M6 e 3 kik 2 k 3 kf 



(3.10) 



To relate this to the familiar result in Chen et al. [16] one uses the fact that E = eM 2 H 2 /c 2 , 
leading to the following expression: 

q #4 



(CfciCfc 2 Cfc 3 )BD 



M*e 2 c 2 s 



1 



1 

k\k 2 k 3 kS 



(3.11) 



As is well-known, the momentum dependence of this contribution to the three-point func- 
tion is well-approximated by the equilateral template, i.e. its magnitude maximizes in 
equilateral triangles. The full three-point function obviously contains additional contribu- 
tions from the other terms in the interaction Hamiltonian. As mentioned earlier, in DBI 
models the above contribution vanishes exactly and another leading order equilateral term 
dominates (in the limit of slow-roll and small sound speed). All other contributions turn 
out to be higher order in the slow-roll expansion [16]. 
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3.3 Bispectrum corrections due to a modified vacuum state 

Next consider the following Bogolyubov transformation of the initial state, 



Uk -> a k v k + P k v* k , 



(3.12) 



with \a k \ 2 — \P k \ 2 = 1. We will be interested in small departures from the usual BD 
vacuum, so |/3fc|/|afe| <C 1, as required by power spectrum constraints. For the BD vacuum 
(a k = 1), our results will of course reduce to (3.10). Now let us suppose (3 k ^ and 
compute the corresponding bispectrum to lowest order in (3 k . This will result in one of 
the modes effectively having negative energy and will carry along a corresponding factor 
of f3 k . As emphasized previously, we will have to introduce an initial time cutoff since we 
expect particle interaction before the end of inflation to dominate the bispectrum, which 
could become arbitrarily large if the earliest modes have infinite time to grow (before they 
cross the horizon). To preserve scale-invariance, we introduce a fixed physical cutoff scale 
and determine the intial time as kc s rjo = A c /H* . Of course, scale-invariance also requires 
that the complex Bogolyubov parameter depends at most only weakly upon the comoving 
momentum k. The starting point will therefore be a constant complex parameter (3 = 
\P\ exp (z<5) that parameterizes the departure from the BD state. Under these assumptions, 
the correction to the three-point correlation function eq.(3.11) is 

H G 



(CfciCfc 2 Cfc 3 ) 



nBD 



Tie 



(2vr) 3 ^E 



kj 



S I 



1 



+ 2A 



-cl x 



drj' 1 
a 2 k\k2k^ 



^ikj C s Tj 



(3.13) 



where kj = kt — 2k j. The integral is evaluated to be equal to 
r-o 



tt2 /-0 _ tt2 

£ / drf r, ,2 e ik ^' 



kik 2 k 3 



kik 2 k 3 ^ 
j 



2 e ikjC s rio 

h i- 

k j C s 



kjC s 



2 , • 2r 7o 

KjC s 



U2 2 
KjC S) 



(3.14) 

The first term is similar to our previous result, with the substitution of kj <-> kt- The 
remaining terms are a result of the presence of a cutoff. The final result is 



3/nBD(l) 



1, 



e 2 cf 



E 



clrj 2 cos(kjc s r] + 5) 
2 jfe. 



c s r]o 



^ kik 2 k 3 
sm(kjC s r]o + 5) cos 5 — cos(kjC s r]Q + 5) 



(3.15) 



+ 



Let us pause for a moment to interpret this expression. First, we notice that all terms 
contain a phase shift 5 as a result of the phase in the Bogolyubov parameter. As will become 
clear this phase significantly influences the projection onto the observational templates and 
therefore the constraints on the vacuum state modification that can be derived. 

Secondly, we notice that in the limit kj — > (and unlike previous results [8]) all 
terms scale as 1/kj, which might cause some concern when considering the limit kj — > 0. 
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We certainly expect the result to be finite for any phase 5 and any combination of the 
three comoving momenta and fortunately explicit computation of the limit kj — > in 
the appendix confirms this expectation. The first maximum occurs for kj oc l/c s i]o and 
since kc s r]Q 3> 1, where k would be the largest momentum scale under consideration, this 
effectively corresponds to the enfolded limit kj/k ~ 0. Also note that this localized form 
of enhancement would scale as (fec s ?7o) 3 , which is expected to be very large. Projection 
onto observational templates generically results in the loss of at most one factor of kc s rjo, 
which would then provide us with a quadratic enhancement in the large number kc s r]o- A 
detailed analysis of the projections and constraints will be discussed in section 4. 

After this rather explicit introduction and discussion of the corrected three-point func- 
tion due to the first term in the interaction Hamiltonian we would now like to present 
the complete (leading order) result for the corrected three-point function, for which the 
detailed calculations are presented in the appendix. Using the expression for the power 
spectrum of general single-field inflation 



8ir 2 c. s e ' 



(3.16) 



we can write the leading order behavior, to first order in (3, as follows 



(C kl C k2 Ck 3 ) = (2vr) 7 |/3]P f 2 <5(^fc i )p^ 



2-i3 Z-a3 h*3 



+ £( B W +B <3) _ B (4A 2£ B (3)V 



klCsVO 



-2s + 1 - C^(5) + e_ / ( 5 ) R (6) 



+ .. 



(3.17) 



where the ... represent terms further suppressed slow-roll and in the cutoff r/o- Here the 
different functions B equal 



B 
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1 



1 



■j+l Kj+2 



So 
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Si 



k\ 



(kjC s r]o) 2 kjC s r] 
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+k 2 k 2 +2 (kj - 2k j+1 - k j+2 ) + kj +1 kj +2 (2kj - k j+1 - k j+2 )+ 



(kjC s rj ) 2 kjC s r] 



cos 5 S± 



S 2 



) 



(3.21) 



(kjkj+i + kjkj+2 - k j+1 k j+2 )Si, 
^3) — (k\k\ + k 2 k 2 + k 2 k 2 ) J S\ 



(3.22) 



(3.23) 



and 



S^kjCgrjOiS) 



cos(kjC s r)o + 6) 



S 2 (kjC s r]o,5) 



kjC s r] Q 
sin(fcjC s ?7o + <5) 



kjC s rj 



We have organized the bispectrum into pieces with different (decreasing) powers of kic s i]o, 
i.e. the ratio between the Hubble radius and the largest physical momentum scale at the 
cutoff time 770, i-e. k\c s rio = (k\ / ao) / (H / c s ) which can be as large as 10 3 . For our purposes 
only the first two terms will be relevant as these are both enhanced in k\c s 770. From here on 
we will refer to these terms as the leading (proportional to (kic s r]o) 3 ) and the subleading 
(proportional to (kic s rjo) 2 ) term, respectively. 

On first impression it seems that the momentum-dependence, or the shape, of the differ- 
ent terms in the above bispectrum strongly deviate from the constrained local, equilateral 
and orthogonal templates. Most importantly, all terms rapidly oscillate. Furthermore, 
as mentioned already, the leading-order term maximizes in the collinear limit. In [8] we 
proposed a template for collinear triangles, which unfortunately is only a marginal im- 
provement as compared to the local and equilateral template. Ideally, one would like to 
introduce a template that can (also) incorporate oscillations with an a priori undetermined 
frequency, introducing an additional parameter. It remains to be seen whether oscillations 
in the momentum dependence of the bispectrum are observable at all, as the finite number 
of pixels in the data might make determination of this additional frequency parameter very 
difficult. We will not be concerned with this problem here since we will only be interested 
in projecting onto bispectrum templates that have already been constrained. We hope to 
report on the opportunities for the detection of oscillations in the primordial bispectrum 
in future work. 

4. Observational constraints 
4.1 Preliminaries 

In order to use our results to place constraints on \/3\ we need to relate them to existing 
templates. This is done by computing the dot product between our bispectrum result 
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Figure 2: Here we plotted the 'density' of the different non-Gaussian shapes. As can be seen from 
the color (shading), large density corresponds to lighter color (shading). It is clear from inspection 
that the equilateral template has zero density on the collinear line (spanned by xi + X3 = 1), 
while the leading order term from initial state modifications maximizes in the vicinity of this line, 
resulting in virtually zero overlap (which is already small due to the oscillatory nature of the leading 
term) . 



and the local, equilateral and orthogonal templates, which are standardized and whose 
magnitude has been constrained in the literature [1, 19]. To facilitate this, let us first 
rewrite our result in terms of the dimensionless variables X2 = ^2/^1 and X3 = k^/ki, 
scaling out a factor of /c^ 6 , indicating that the bispectrum is scale-invariant, where the 
comoving momentum magnitude k\ is identified as the largest momentum. Because of the 
triangle constraint the shape function then only depends on a finite (triangle) domain in 



x 2 and X3 space (0 < X2 < 1 and 1 — X2 < X3 < 1). Using these dimensionless variables, 
the leading shape function is then expressed as 



^ieading(Z2,S 3 ,p,<5) = 



X\X2 



1 cos(p(xi + X2 — 1) + 5) sin(p(xi + X2 — 1) + 5) 

2 p(xi + x 2 -l) p 2 (xi + x 2 -l) 2 

(4.1) 



cos 5 — cos(p(xi + X2 — 1) + 5) 
+ 37 : tto + eye. perm. 

p 6 {Xi + X 2 - l) d 



where p = kic s r]o is the enhancement factor we encountered several times before. Having 
identified the shape of the leading contribution we automatically also fixed a particular 
expression for the non-Gaussian amplitude /nl- It is important to note that the above 
shape function implies that a factor of p 3 has been incorporated in the non-Gaussian 
amplitude /nl- 

There are three different observational templates that we will be projecting our signal 
onto. These are the squeezed or local, equilateral and orthogonal templates. In terms of x 2 
and xs, their shape functions are expressed as follows: 

iWi,* 2 ,x 3 ) = 2 + h + i) > ( 4 - 2 ) 



F eq (l,x 2 ,x 3 ) = 6 ( - ^— + _J_ -1 -1 ?_ + perm. ) , (4.3) 



Fort(l, X 2 , S 3 ) =3(^-^2). (4.4) 



X <2p£> ^ 



Each of these templates has a different physical origin. The squeezed template was in- 
troduced by Spergel and Komatsu [36] to parameterize non-Gaussianities resulting from 
a quadratic correction in real space (i.e. a local term), corresponding to the leading non- 
linear coupling between the inflaton field and the curvature perturbation at superhorizon 
scales. The equilateral template was introduced by Creminelli et al. [22] to parameterize 
nonlinear higher-derivative corrections in the inflaton field, whose nonlinear magnitude is 
fixed at horizon crossing. As mentioned, the orthogonal template was recently introduced 
by Senatore et al. [19] to allow for an optimal comparison with exotic non-canonical single- 
field scenarios that are, rather surprisingly, not of the equilateral type. A closely related 
template was also put forward in earlier work [8]. 

In order to make comparisons between these different shapes, in [22, 23] dot products, 
cosines (the normalized dot product) and projection factors (to be introduced later) were 
defined to be 

k 4 k 4 k A 

dkidk 2 dk 3 12 3 F x (fa, k 2 , fa)F Y (ki, k 2 , k 3 ), 



F X ^F Y = 



k. 



COS(i?X ' Fy) 55 {F x *Fx)l^y*F Y y/r ^ 

In the case where both Fy and Fx are scale-invariant we can remove the integral over one 
variable, (i.e. ki). In a previous definition of the dot (or star) product [22] the factor 1/kt 
was omitted. It was argued in [23] that the addition of this extra 'weight' reduces the differ- 
ence between the three-dimensional dot product for the bispectra and the two-dimensional 
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dot product for the respective multipole equivalents. Since one actually measures the three- 
point function in two-dimensional multipole space, it is the two-dimensional dot product 
that one is interested in computing to derive constraints. Since the two-dimensional mul- 
tipole bispectrum and its corresponding dot product are computationally challenging (and 
analytically ill- understood) , it would indeed be a great improvement if a three-dimensional 
dot product could be introduced that guarantees reproduction of the two-dimensional dot 
product in multipole space. This is what the above dot product (including the 1/kt) is 
supposed to do, removing the need for a fully fledged two-dimensional analysis. 

We will take this opportunity to clarify why the assumption of scale-invariance pro- 
duces a difference with the old dot product which is (approximately) only a factor of 
1/2. Note that the weight function 1/kt is given by 1/(1 + 22 + 23) (after scaling out 
k\), and the integration limits are given by < 22 < 1 and 1 — 22 < 23 < 1. Conse- 
quently 1 + 22 < I + 22 + 23 < 2 + X2 or 1 < 1 + X2 + 23 < 3, resulting on average in 
kt = 1 + 22 + 23 ~ 2, explaining the factor 1/2 difference. When computing the projec- 
tion factor, a ratio of dot products that we are ultimately interested in, this factor 1/2 
drops out and one ends up with exactly the same result using the former definition of the 
projection factor. We should emphasize that this is only true under the assumption of 
scale-invariance. Once this assumption is dropped one cannot scale out the k\ explicitly 
and one has to perform the full three-dimensional integration in A;-space. Since the dot 
product without the weight 1/kt can be computed analytically we will actually be working 
with the standard definition of the three-dimensional dot product (see appendix eq. (B.l)) 
and assume the results will (approximately) equal the two-dimensional multipole space dot 
product. 

Since we expect our bispectra to have poor overlap with the templates, and the dot 
product between the theoretically predicted bispectra is very hard to compute analytically, 
we do not consider the cosine between the spectra. However, to get an indication how well 
(or how poor) a template does as compared to each other, we can define the ratio between 
two cosines 



K(X,Y) 



cos (F x , F) 



cos(Fy,F) 



Fx-F 



Fy ■ F 



where F represents the theoretically predicted bispectrum, and Fx and Fy represent dif- 
ferent bispectrum templates. If TZ(X, Y) > 1, Fx is a better template and if TZ(X, Y) < 1 
the opposite is true. 

Figure 2 shows the density of each shape which gives a qualitative indication of where 
the shapes peak as a function of the comoving momenta. We can immediately see that the 
orthogonal template naively appears to 'best' match the qualitative features in our signal 
(at least for 5 = 0), an intuition that will be (partially) supported by calculations of the 
dot product. As an instructive example, let us consider the cross-correlations with the 
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leading shape function eq.(4.1), which can be analytically determined, 

w v - 3cosS -i- sin6 ±n(rr^ 

r local • -Heading — T 1 9 r U(p ), 

4 p p z 

sin S 

Feci ■ -Heading = 2 s h C(y~ 3 ), 

-Tort • -Heading — ~ r O s h U(J> J. 

2 p p z 

Consequently both the dot product with the local as well as the orthogonal template peak 
for 5 = 0, while the dot product with the equilateral template is ir/2 out of phase. Note 
that we extracted three powers of p into the amplitude, so the dot products should be 
multiplied with p 3 to determine the strength of the enhancement that remains. The reason 
for this reduction in the level of enhancement is simply that the original enhancement was 
localized in the (leading) three-point function and integrating over the full domain will 
smear out this enhancement, typically reducing it with one power of p. Clearly, the dot 
product with the equilateral template is additionally suppressed in p as compared to the 
local and orthogonal template (for 6 = in fact all enhancement is lost). Visual inspection 
of the leading-order term in the density profile of figure 2 indicates that the small dot 
product with the equilateral template is explained by the fact that this term peaks in the 
collinear limit, where the equilateral template actually vanishes. As mentioned before, 
this is precisely where the orthogonal template peaks, and therefore it is not a surprise 
that it would have a larger dot product. A significant feature of the orthogonal template 
is that it is not positive-definite within the triangle domain. As such it shares a feature 
with the oscillating contributions, possibly resulting in an even larger dot product. For 
comparison we find TZ(Fi oca x, F ort ) ~ 0.42 for 5 = which quantitatively confirms the 
orthogonal template has the best overlap with the theoretical spectrum 1 . 

We conclude that for generic values of the phase 5 (5 ^ vr/2), the local and orthogo- 
nal templates would provide the best constraints if the shape function eq.(4.1) would be 
a good approximation to the full three-point function. As one should expect, the con- 
straints resulting from these expressions will also depend on the tightness of the existing 
observational constraints on the corresponding /nl, which will favor the local template. 

For completeness we should also consider, in addition to the term (4.1) which is a 
unique consequence of introducing a vacuum state modification, the non-Gaussian signature 
as a consequence of non-standard c s (i.e. O(/3 )). As computed in [14], the leading terms 
are given by 

1 f I 2A\ Skfkjkj 



where we neglected the corrections due to time-variation of various parameters as well as 
corrections to the solutions of the field equations of motion. Both of these contributions 

1 A& a reference, the overlap between the local and equilateral shape is cos(i*i oca i, F cq ) ~ 0.41. 
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have a /c-dependence that strongly resembles the equilateral template, 



cos(F cq , B x ) ~ cos(F cq , B c ) « 0.95. (4.9) 

This contribution is therefore expected to add a term to the equilateral non-Gaussian 
amplitude, in addition to the modified vacuum state effect. It will of course be interesting 
to see how these contributions compare. 

4.2 General inflationary models with a small speed of sound 

We can use the dot product, i.e. the correlator, to define projections, which measures the 
relative "leakage" of one shape function into one of the available observational templates 

A(F X ,F Y )^^-^. (4.10) 

Using the normalizations (self-dot products) of the different templates that are given in the 
appendix, it is straightforward to compute the projection factors for the different templates 
with the leading shape behavior of the corrected three-point function due to a modified 
vacuum state. For large p this is well-described by the shape function that was discussed 
in the previous paragraph. The results are 



A{ f leading; -f local) 



1 



176.5 



3 cos 5 sin 5 

4 p p z 



(4.11) 



A(F Ieading ,F eq ) « A™*, (4.12) 



A (f leading; -fort) 



3 



13.8 



1 cos 5 sin 5 
+ 



(4.13) 



6 p p 2 

A(B x ,F eq ) « 0.01, (4.14) 
A(B c ,F cq ) « -0.05. (4.15) 

From the corresponding amplitude of the three-point function, including the p 3 enhance- 
ment factor, one then derives the following expressions for the local, equilateral (including 
contributions due to a small speed of sound) and orthogonal non-Gaussian amplitudes: 

- ^ (HI cos* + Ip-^S) + l/V, (4-16) 

/nl = -3 (^f l/3bsin5 + 0.04^ + ^l^sin^-O.Ol) , (4.17) 

^ ~ its (1™ sS + 5 p- 1s ™ s ) - 1 + Wp 2 - ( 4 - 18 ) 



Again, in deriving these results we have assumed, based on the arguments presented earlier, 
that the three-dimensional dot products are good approximations of the two-dimensional 
dot products of the multipole bispectra. Also note that the two terms in the equilateral 
non-Gaussian amplitude that are not proportional to \/3\ are a direct consequence of a 
small speed of sound, which produces an enhanced equilateral non-Gaussian component 
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which cannot be neglected for 5 = 0. The local and orthogonal contributions are, however, 
uniquely due to the modified vacuum state. We are now in a position to use these predic- 
tions to put constraints on the parameter space, consisting of c s , \/3\, 5 and p. The existing 
constraints on f l °™ 1 [37], /^ uil [1,19] and [19] are given by 

-4 < /i?£ al < 80, 



'NL 

-125 < /* q L uil < 435, 
-369 </gf< 71. 



Based on these results let us now consider two extremes: 5 = (real (3) and 5 = tt/2 
(imaginary (3). 

For 5 = 0, eq. (4.12) vanishes, which implies the linearly enhanced equilateral con- 
tributions proportional to \(3\ vanish. In that case we can actually use the equilateral 
constraints to provide a bound on the speed of sound, similar to that found by Senatore et 
al. [19], but as we will see this is not even required to derive an interesting constraint on 
\/3\. One should also realize that the 5 = case describes the generic order of magnitude 
expectation, in the sense that as long as 5 ^ tt/2 the local and orthogonal contributions 
remain p 2 enhanced, which is the most important property as far as the constraints are 
concerned. The constraints from the orthogonal template and the local template are very 
similar. Since \P\ is positive-definite this implies only the positive constraint remains (the 
negative constraint can be used for 5 = tt). Note that the projection factor between the 
local and the leading shape function has a negative sign, so the final constraint comes from 
the (negative) lower bound on local non-Gaussianities /]^ al , which is given by 

WS^d-l+^'.W-O). (4.19, 

To make this more concrete, let us assume that the ratio 2A/S <C 1/c 2 in the small 
sound speed limit. The constraint then reads \/3\ < 6 x 10 2 c 2 /p 2 and after realizing that 
p = kic s riQ = A c c s /H this implies the interesting, sound speed-independent, constraint 
\(3\ < 6 x 10 2 H 2 1 A 2 . This could easily give rise to a constraint on \(3\ at the level of 10 -4 , 
independent of any specific model for \/3\. In particular it seems to exclude Bogolyubov 
parameters linear in H/A c [31] 2 . Therefore in the small speed of sound limit (and 2A/S <C 
1/c 2 ), we obtain an impressive constraint on \/3\. As long as the phase 5 ^ 7r/2 this 
conclusion is unaffected. For values 5 > tt the sign in the leakage factors is altered and the 
best constraint is obtained from orthogonal non-Gaussianities. Due to the larger number 
bounds on this /nl> this results in a constraint that is one order of magnitude below the 
best constraint. 

Now let us briefly consider the special case 5 = tt/2. For the local and orthogonal 
contributions this means the leading p behavior is one order of magnitude less, which 



2 Strictly speaking, starting with the assumption that the Bobolyubov parameter is proportional to 
which introduces another factor of c s , one would obtain a lower bound on the ratio of depending on the 
speed of sound which turns out to be in conflict with constraints from the power spectrum in the limit of 
a small sound speed. 
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Figure 3: Constraining contours for different values of p = c s kir]o and two values of 8. The 
constraints for intermediate values of 8 should lie somewhere in between. For values tt/2 < 8 < 3tt/2 
one has to consider the opposite (in sign) constraints on the various /nl- The constraint for p = 1000 
and 8 = is too tight to be visible in this graph. For 8 ~ tt/2 and p — 1000 the constraint is still 
weaker than the constraint for (5 = and p = 100. 



implies that the equilateral contribution is now of the same order. Looking at the equilateral 
contribution we conclude that we can neglect the equilateral contributions from the zeroth 
order bispectrum (i.e. without initial state modification) as long as \(3\ > 0.1 /p. If the 
Bogolyubov parameter is smaller than that we can actually use the equilateral contribution 
to put a constraint on c s . This time, the strongest constraint is derived from the orthogonal 
template and reads 

,,,<^(i- 1 + |)-\<^/2). (-) 

Obviously, this constraint is much weaker and can be written as \j3\ < 2 x 10 2 Hc s / A c , in the 
limit of a small speed of sound and 2A/E <C 1/c^. Since the enhancement in p is reduced in 
this case one should worry about subleading contributions since they might be comparable 
in magnitude. Upon checking this one finds that for 5 = tt/2 one should indeed take into 
account both the subleading and the leading order contribution. In the next section we 
will derive the projection factors for the subleading terms, since these are also the relevant 
terms in DBI inflation for which (l/c^ — 1 + 2S/A) vanishes identically. 
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4.3 DBI models of inflation 



For DBI models (IR or UV) the leading-order term of the vacuum state corrected three- 
point function vanishes. Consequently we will need to compute the dot product with the 
next-to-leading order terms. The next-to-leading contribution goes as (from eq.(4.21) and 
eq.(4.21)) 

where we have 2 distinguishable contributions, 
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As was the case for the leading contribution, both terms are finite in the collinear limit (as 
shown explicitly in the appendix). Both contributions are of the same order of magnitude 
over the full triangle domain, so both have to be included in our analysis. One might argue 
that the cosine in both of these functions is suppressed and can be neglected. However, for 
an analytical treatment of the projection factors it turns out to be convenient to work with 
the full expression because of its finite collinear limit. Moreover, as we will see shortly, the 
phase 5 significantly influences the final results for the projection factors. 

As before, let us rewrite these contributions in terms of the dimensionless variables %i 
and x 3 . We will identify the shape of the first sub-leading contribution as follows 



3 1 



4 X2X3 



1 \ sm(p(x 2 + x 3 -l)+5) . 421 , 

X2 X 3 ) p(x 2 + X 3 -l) 



1 - 2x 2 + x 3 1 + x 2 — 2x 3 \ cos 5 cos(p(x2 + X3 — I) + 5) 
x 2 + x 3 - 2 -\ 3 h 



sl(l) 

( 1-2 

_l_ I <r-„ _l_ ry 9 J 

x\ x\ ) P 2 {X2 + X3 — 1) 

plus cyclic permutations. The shape of the second sub-leading term equals 
1 1 



+ (4 ' 22) 

cos 5 — cos(p(x2 + X3 — 1) + 5) 



((1 + x 2 + x 3 )(x 2 x 3 + (x 2 + x 3 - l)(x 2 + x 3 - X 2 X 3 )) 



p 2 (x2 + x 3 - I) 2 



plus cyclic permutations. These identifications of the shape function imply we have ex- 
tracted a factor of p 2 into the non-Gaussian amplitude. 
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Figure 4: The shape of the first subleading contribution (eq.(4.22)) for different values of 6. For 
S = the shape of this function is very similar to the leading order term. Once 6^0 the shape 
picks up a significant local feature which leads us to expect a significant overlap with the local 
non-Gaussian template. Since the oscillations arc still present, we also expect a comparable overlap 
with the orthonormal template. As for the second subleading contribution (eq. (4.23), not shown) 
we find similar behavior, but opposite in sign, resulting in a small but non-negligible cancellation. 



Before computing the dot products let us have a crude first look at the shape functions. 
They appear quite similar to the leading order shape, but somewhat surprisingly, as soon 
as 5 deviates from zero the sub-leading contributions change qualitatively. On the left in 
figure (4) we have plotted eq. (4.22) for p = 100 and 5 = 0, corresponding to a rapidly- 
oscillating function that maximizes in the collinear limit. On the right however, we plotted 
the same shape function with the phase S = 3ir/2, which suddenly appears to resemble 
the local shape. Interestingly, the appearance of local type features was also noted in the 
different context of multi-field DBI models (assuming the standard Bunch-Davies vacuum 
state) [38]. Based on this observation one would expect the dot product with the local 
shape to increase as a function of 6 (as opposed to what happens for the leading shape 
function). Of course, the collinear feature is still present, but simply swamped by a strong 
local feature that appears for 5 = 3n/2. In the previous analysis of the leading term we also 
found that the dot product with the orthogonal template was ir out of phase as compared 
to the dot product with the local template and that it was comparable in magnitude. We 
would expect to find something similar in this sub-leading DBI case. 

Because of the clear presence of local non-Gaussian features for 5 > 0, which are 
actually divergent in the X2 = or X3 = limit, we have to introduce a cutoff when 
computing dot products, the same cutoff used to normalize the local template. This cutoff 
is necessary to reflect that we cannot observe infinite wavenumbers in the CMB. It is 
given by fcmin/^max ~ 10 -3 , the largest observable physical scale divided by the smallest 
observable physical scale. So we will require that X2,x% > e, with the cutoff e equal to 
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Figure 5: The dot product between the subleading (F s \ = f s i(i) + Fsl<2)) an d local template 
(multiplied with p) for two values of p. Although the dot product is much smaller for S = it 
is non-zero for both p. However, there exist values of the phase for which the dot product (and 
consequently the leakage or fudge factor) is exactly zero. For such S the local template would not 
be suitable (or at least limited) for detecting any deviation from a BD state (in a DBI model of 
inflation) and we will have to use the much less capable equilateral template. 

e = 10 -3 . 

We find that there is substantial overlap with both the local and orthogonal templates. 
For the dot product with the local template we obtain 

+ Fsl(2)) ■ ^locai = ^ (9 + 8 log e) sin 5 + [2 cos(5 + 2pe) - 2(1 - Trpe) cos 5 

+ 4pe(ci(2pe)sin5 + si(2pe)cos<5)] +0{p~ 2 ,e). (4.23) 

Here ci = cosintegral and si = sinintegral. Note that since e ~ 10~ 3 we should be care- 
ful identifying the second term as p~ 2 suppressed. For small values of p the first term 
dominates. First of all, remember that a factor of p 2 was scaled into the non-Gaussian 
amplitude /nl, explaining the appearance of the leading 1/p suppression. Considering the 
case S = 7r/2, e = 1CT 3 , we find that (i^m + F s i(2)) " ^local ~ — 31/p. Properly comparing 
this to the dot product of the local template with the leading order shape function (which 
requires scaling in one factor of p, to match the normalizations of the amplitude), we find 
that that the overlap between the subleading term and the local shape is almost 42 times 
better, which is rather surprising. 

In figure 5 we have plotted the dot product for two values of p as a function of 5. It 
clearly shows that even though the dot product is smaller for 5 = it is still relatively large 
compared to the dot product between the leading term and the local template. There exist 
values for 5 for which the dot product is exactly zero. For those values of the phase the 
local template would not be useful to derive constraints on vacuum state modifications. 
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Figure 6: The dot product between the subleading and the orthogonal template (multiplied with 
p) for two values of p. This dot product is almost completely dominated by the sines and therefore 
almost zero for (5 = 



We were unable to obtain an analytical expression for the dot product with the equilat- 
eral template and the subleading bispectrum. Instead we performed a numerical analysis, 
which is discussed below. Next, we compute the dot product with the equilateral and 
orthogonal templates. Semi-analytically we obtain from best fitting: 

sin S loe v 

GPrf(i) + Fsl(2)) • ^ort ~ 5.3— + 11-3-4^ cos(5 - 0.2 - 2Ay/p). (4.24) 

We compute that H(Fi oca \, F OTt ) ~ 1.71 for 5 = tt/2, indicating that, as expected, the 
local shape has the best resemblance with the subleading bispectrum. 

As in the previous section we use these to obtain the projection factors and then 
translate these into expressions for the various /nl (which introduces a factor of p 2 that 
was absorbed in the non-Gaussian amplitude): 



r local ^ 1 5 
JNL 176.5 6 



10 3 

(9 - 32 log 10) sin 5 + (2 cos(5 + 2 x 10 3 p) - 2(1 + vr x 10" 3 p) cos 8 

P 

+ 4 x 10~ 3 p(ci(2 x 10~ 3 p) sin 5 + si(2 x 10~ 3 p) cos 6))] - lj \/3\p (4.25) 

and 

/nl - ^3 ( 5.3 sin 5+ 11.3^ cos(<5- 0.2 -2.4^)1 I ~ - 1 1 \/3\p. (4.26) 

Once more we distinguish two values of 5. For 5 = tt/2 and in the small sound speed limit, 
using that p = A c c s /H, the best constraint comes from local non-Gaussinaties and reads 

14/7 r 

\P\ < -^A (S = */2), (4-27) 
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which is substantially tighter than the constraint from leading order contributions for the 
same value of the phase 5 (eq. (4.20)). For all the values in between (i.e. < 5 < 7r/2) the 
constraints will in general be somewhat weaker. Once the sign of the leakage factor changes, 
the constraints loose about an order of magnitude, which still exceeds the constraints from 
the leading term. From eq. (4.25) we find the following constraint for 5 = 

l/?l<^, (5 = 0). (4.28) 

For 5 ~ 2.9 eq. (4.27) is zero and therefore not useful to place any constraints on 
Similar to the dot products with the leading order terms the dot product with the 
equilateral template is approximately ir/2 out of phase with the orthogonal and local ones. 
Obtaining an analytical solution for the dot product between the equilateral template 
and the subleading non-Gaussianities turned out to be technically involved, instead we 
computed these numerically for different values of p and 5. These are shown in figure 7. 
What we find is similar to what we found for the leading order term: it is indeed out of 
phase with the other two dot products, nor is it enhanced in p (only log-scale). These 
observations imply that although in principle one could use the equilateral template to 
constrain \@\ in the vicinity of those values for which the other two templates would not be 
able to place a constraint at all, due to the poor overlap with equilateral non-Gaussianities 
any constraint on \(5\ would be very weak. It must be noted that the range of values for 
5 for which the constraints are weak is limited and it appears, due to small corrections, 
that the zero points of the the dot products between local and subleading and orthogonal 
and subleading do not exactly coincide. As a result, one might still be able to obtain a 
slightly better results using the orthogonal constraints. Since the equilateral template is 
only log-scale enhanced by the cutoff, given the constraint on |/3| from the power spectrum 
(roughly \(3\ < 10 -2 ), the signal in the equilateral template will sometimes be dominated 
by the 'regular' DBI contribution, which would allow us to put a constraint on c s thereby 
breaking the degeneracy between c s and \/3\. Unfortunately, this is only true for a small 
range in 5, effectively only for those values of 5 for which the equilateral contribution from 
vacuum modification is zero in any case. From best fitting our numerical results, we have 
been able to deduce that 

(F sl{1) + F sl{2) ) • F cquil ~ -10.8cos(<5 - 0.3) x (4.29) 

Consequently, can be written as 

/nI" 1 - -|(0.05 + 10.8cos(,5- 0.3) x \P\logp) (± ~ l) • (4-30) 

For DBI models of inflation the speed of sound c^ 1 can easily be O(10). In such a 
scenario for generic values of the phase 5 the constraints become 

101 £ (S = n/2), (4.31) 
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Figure 7: Numerical results for the dot product (multiplied by p) between the equilateral template 
and the subleading terms. The dot product gets smaller for larger values of p, indicating that the 
constraints using this template are not improved for larger values of p (recall that one p was absorbed 
into the definition of the associated /nl)- In addition, the dot product is already small to begin 
with, indicating that the equilateral template might be useful to constrain DBI inflation itself, but 
is rather poor in constraining DBI inflation with modified initial states. Note that for p = 500 
the numerical computation is at its limit; for p values beyond this the numerical precision is not 
sufficient to reproduce the expected (cosine) oscillations. This holds especially for values that are 
close to zero. There was no need for a cutoff (e) in these calculations as the dot product was not 
divergent (as expected) for small values of X2 and X3. 



which implies that for realistic values of the cutoff \[5\ < 1CP 3 . This bests the constraints 
from the power spectrum by one order of magnitude. As noted however, there exists 
a limited range in 5 for which the constraints are weaker due to the poor overlap with 
existing templates. Unlike the constraint coming from the bispectrum, the constraints on 
\(3\ obtained from the power spectrum do not depend on the cutoff scale or the phase of 
the Bogolyubov parameter. As such, the bound on \(5\ for these values of the phase is best 
constraint by the power spectrum. 

5. Conclusions 

We have analyzed three-point correlators in single- field inflation with a small speed of sound 
after introducing a small departure from the standard Bunch-Davies vacuum. Our aim was 
to derive the best possible constraints on slightly modified vacua in the specific context 



and 




(4.32) 
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Figure 8: Constraints on |/3| as function of 1/cJ — 1 derived from the dot product and the 3a 
levels of local non-Gaussianities for different values of p and S. These constraints are better than 
the constraints coming from measurements of the power spectrum. However, this requires several 
assumptions due to the appearance of 4 independent variables (p, c s , \J3\ and S). p could be measured 
independently if one could be able to measure the frequency of the oscillations in the bispectrum. 
c s could be determined if one could deduce what type of DBI model initiated inflation. As such, 
one could break down some of the degeneracy. 

of small sound speed models like DBI inflation. In previous work [8] we concluded that 
the bispectrum is an excellent probe of vacuum state modifications due to enhancement 
of the non-Gaussian signal by a factor proportional to (powers of) the ratio of the cutoff 
scale and the Hubble parameter p = A c /Hc s , which could easily be as large as 10 3 . We 
have shown that in the small sound speed limit the powers of p that appear can be as 
large as p 3 , which adds one factor of p as compared to the c s = 1 model. In addition 
the level of non-Gaussianity in small sound speed models is already large to begin with. 
Therefore it is not surprising that the constraints that can be derived in such models are 
in general stronger. The computation of the contributions to the local, equilateral and 
orthogonal template confirmed this expectation. In general single field models with a small 
speed of sound the localized enhancement is proportional to p 3 , and after projecting onto 
the local and orthogonal template one factor of p is lost, ending up with p 2 enhancement. 
The inclusion of an arbitrary phase in our analysis in the Bogolyubov parameter turns out 
to be able to reduce the enhancement with one factor of p for very specific values of the 
phase. For general single field models with a small speed of sound this implies that for 
generic phases the constraint on the absolute value of the Bogolyubov parameter equals 
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\P\ ^ 6 x 10 2 H 2 1 A 2 , which is actually independent of the speed of sound. This strong 
constraint is in fact excluding Bogolyubov parameters linear in H/A c in the small sound 
speed limit. To be precise, assuming \/3\ oc ^- one obtains a (c s dependent) lower bound on 
the ratio which in the limit of a small sound speed would violate the power spectrum 
constraint. For the specific value of the phase 5 = tt/2, the subleading terms are of the 
same order and have to be included in the analysis. Similarly, it is the subleading terms 
that govern the constraint for DBI models, for which the leading term identically vanishes. 

In the context of DBI models we performed the same analysis, with a reduced en- 
hancement of p 2 in general. After projecting onto the different templates and enforcing the 
observational bounds this resulted in a constraint on the absolute value of the Bogolyubov 
parameter \(3\ < 14H/(A c c s ), which is linear in H/A c and explicitly depends on the speed 
of sound. This subleading contribution is the one that should also be taken into account 
for the specific value of the phase 5 = tt/2 in the general single field case, since this bound 
is much stronger than the bound from the second order terms in the leading expression. 
For general single field inflation with a small speed of sound this results in a very tight 
bound on allowed values of \(5\ ranging from 10~ 5 to 10~ 3 . For DBI models of inflation 
the constraints are only determined by the subleading term as the leading term vanishes. 
Again, the inclusion of the phase of the Bogolyubov parameter has the effect that for the 
specific value of 5 = the constraint weakens slightly, equal to \[3\ < 72H/A c c s , but this 
time remains at the same order in p. Overall we conclude that for generic values of 8 and 
realistic values of the cutoff this still implies \/3\ < 1CT 3 — 1CT 2 . For very specific values of 
the phase, the bounds could vanish altogether. For these special phases one would have to 
rely on the power spectrum to set a constraint on \(5\. 

Our main message is that for single field models with a small speed of sound, including 
DBI, the constraints on vacuum state modifications from the bispectrum for generic values 
of the phase are already better than the power spectrum constraints, relying on the available 
templates. The main reason for this is the huge (localized) enhancement that appears 
after introducing a modified initial state and a corresponding cutoff that regularizes the 
relevant integrals for the bispectrum calculation. Improved templates, which would have 
to be sensitive to the oscillatory nature of most of the contributions to the bispectrum, 
are guaranteed to tighten these constraints significantly. As was previously noted [8], 
the oscillatory nature of the signal in momentum space suggests that specific, perhaps 
observable, features could appear in three- and two-dimensional position space. In any 
case it would be worthwhile to generalize the range of available non-Gaussian shapes that 
can be compared to the data, including oscillatory signals, on which we hope to report in 
the future. 

Our results confirm that higher- derivative corrections are excellent probes for depar- 
tures from the standard Bunch-Davies vacuum state. Throughout this paper we assumed 
that the combination | A; 1 770 1 is independent of the actual comoving momenta involved and 
equal to A c /H, in the spirit of the New Physics Hypersurface approach to vacuum state 
modifications. The reason for this was scale-invariance of the bispectrum, which we relied 
on to allow for comparison with the available (scale-invariant) template shapes. Fixing 
770 instead, as one would do in a Boundary Effective Field Theory approach to vacuum 
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state modifications, immediately results in a scale-dependent bispectrum. It would be in- 
teresting to study such scale-dependent scenarios and determine to what extent (future) 
analysis of 3d large scale structure or 2d CMB data can constrain bispectrum departures 
from scale-invariance [21]. 

As reported, the bispectrum or three-point function is extremely sensitive to initial 
vacuum state modifications in the presence of a higher-derivative operator, and there is 
no reason to think this could not similarly be true for all higher n-point functions. A 
more general perturbative analysis, including higher n-point functions [39,40], might lead 
to a hierarchy of (theoretical) constraints on vacuum state modifications, perhaps ruling 
out departures from the standard Bunch-Davies state altogether. Indeed, for the four-point 
function (trispectrum) the effect of a vacuum state modification was recently studied in [40] , 
again showing a strong sensitivity. We hope that ongoing future work in this direction can 
help us further understand and identify the phenomenological and theoretical constraints 
the inflationary vacuum state has to satisfy. 
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A. Full Bispectrum 



In this section we will compute the full bispectrum coming from all terms in the inter- 
action Hamiltonian and for general but /c-independent (3 with an arbitrary phase, i.e. 
(3= \0\exp(iS). 

We have already computed the leading-order contribution. In the colinear limit kj — > 
this reduces to 
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Performing a similar analysis as done for the first term on the remaining terms in the 
the interaction Hamiltonian yields for the second term: 
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In the colinear limit kj ->0we obtain 



Hin<C fcl Cfc a C*3>nBD(2) = ^(2vr) 3 ^( e -3 + 3 C 2) ( 5(^fe 



'A 

e 2 ci 



(kf + hk 2 + k\) ((fcifc 2 (fci + k 2 )yoc s cos 5 - 2{k\ + k x k 2 + k\) sin£) 

k\k\{k x + fc 2 ) 



plus cyclic permutations. There are two things of note here. First of all, it is clear that 
the subleading term in this limit has two contributions that are out of phase. Secondly, 
the term proportional to sin 5 is only linearly enhanced. Therefore, for large values of the 
cutoff this term is negligible. Consequently, this implies that the leading and subleading 
are out of phase in their enhancement. Although this observation seems preliminary given 
that we are considering the colinear limit, it turns out as shown in section 5, that this is 
true for the full domain. 
The third term yields 
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This can be integrated to find 
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The contribution from the third line would be infinite in the limit r\ — > and therefore 
would produce a divergent result. However, the limit of the conformal time to zero is not 
an actual limit since this would suggest we would have infinite time to observe all scales. 
This limit is usually taken for convenience since it turns out that in most cases it does not 
add power to the 3-point function. In order to constrain this term one should take the 
limit Tj — ► Tj* where rj* is the conformal time at CMB formation. In this way, this terms is 
finite and negligible compared to all other contributions. We can investigate the behavior 
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of this bispectrum in the limit kj —* to obtain 

Jim (C fcl Cfe,Cfc 3 >nBD(3) = h2it? H ^{e-2s + l-c 2 s )5( i £h) 

(kf + kik 2 + fcf) ((fcifc 2 (fci + k 2 )r) c s cos <5 - 2(kf + fejfo + fc|) sin 5) 

plus cyclic permutations, making this limit exactly similar and opposite in sign to the 
previous collinear limit. In effect it implies that for all values of the phase 5 the collinear 
limit is zero for DBI models of inflation. 

For the last term in the interaction Hamiltonian we obtain 
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B. Template normalization and previous results 

The projection between two shapes Fx and Fy is determined via the dot product 

Fx ■ F Y = / dx 2 dx 3 F x {x 2 ,x 3 )F Y (x 2 ,x 3 )xix3, (B.l) 

first introduced in [22]. This was previously revised and replaced with the improved dot 
product (eq. (4.5)) in [23]. Although this revised 'correlator' improves the results compared 
to the observed 2-dimensional projected dot products, in the case of scale-invariance the 
differences are marginal. The limits of integration are given by the triangle domain, i.e. 

hi = 0) which corresponds to < x 2 < 1 and 1 — x 2 < x 3 < 1. One can then define 
a cosine between two different shapes to 'measure' how similar two shapes with different 
comoving momenta are 

™ fefr)5 (ftWw5' (B ' 2) 

In previous work we calculated the normalization factors F\ oca i • -Fiocal = 176.5 and 
F eq ■ F eq = 7.9. The normalization of the newly proposed orthogonal template F ort ■ F OTt = 
13.8. Note that the normalization of the local template is actually infinite if one would 
integrate over the full triangle domain. This is because most of the signal comes from scales 
k — > which correspond to perturbations with an infinite physical length scale. These are 
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not observable in the CMB and therefore one truncates the limit of integration such that 
only observable scales are included, which we take to be %2 > 0.001, as this is the ratio 
between the smallest and largest observable scales. 
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